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Q>^ I Four different types of free energies are computed by both thermodynamical Bethe 

Ansatz (TBA) techniques and by weak couphng perturbation theory in an integrable 
one-parameter deformation of the 0(4) principal chiral cr-model (with SU(2)xU(l) 
^. symmetry). The model exhibits both 'fermionic' and 'bosonic' type free energies 

r-| ! and in all cases the perturbative and the TBA results are in perfect agreement, 

strongly supporting the correctness of the proposed S matrix. The mass gap is also 
is-H computed in terms of the A parameters of the modified minimal substraction scheme 

j_| ■ and a lattice regularized version of the model. 
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1. Introduction 

We shall study the following one-parameter family of 1+1 dimensional a-models 
described by the Lagrangian 

^0 = -^{LlLl + LlLl + (1 + g)LlLl) , (1.1) 

where 



-L^a" = G-^d^G , G e SU(2) , (1.2) 
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g isa (real) parameter and a" stand for the standard Pauli matrices. The Lagrangian 
( |1.1|) can be interpreted as a deformation of the SU(2)xSU(2) (or 0(4)) symmetric 
nonlinear a-model (NLS) by the parameter g. It possesses an SU(2)lxU(1)r (global) 
symmetry. We shall assume in the following that the deformation parameter, g, 
satisfies 

-l<g<0, (1.3) 

because with g being in this range the model, ( |1.1| ), is asymptotically free (AF) (see 
Appendix D). At the classical level the deformed model, (|1 . 1|) , interpolates between 
the SU(2)xSU(2) {g = 0) and 0(3) {g = -1) NLS models. It has already appeared 
in the Bethe Ansatz approach of Polyakov and Wiegmann to solve the SU(2) x SU(2) 
and 0(3) NLS []T], Q. The deformed model ( |1 . 1| ) is known to possess a Lax-pair and 
is believed to be integrable both at the classical and at the quantum level 0. Its S" 
matrix has been proposed long time ago fl], @] and at present a considerable amount 
of evidence corroborates the belief that the deformed SU(2)xU(l) symmetric a- 
model is quantum integrable indeed 0, ||, §. Its spectrum is thought to contain 
two massive doublets whose scattering is described by the tensor product of an 
SU(2)xU(l) symmetric solution of the bootstrap ^-matrix equations: 

s{e) = s'^'^\e) ^ s'^p\e) , (i.4) 

where S^^\9) is the Sine-Gordon (SG) S-matrix depending on the parameter p. In 
Eq. ( |1.4|) the limit p ^ oo corresponds to an SU(2)xSU(2) symmetric S'-matrix. 

In this paper we carry out a rather thorough consistency check on the proposed 
5- matrix of the one parametric deformation of the 0(4) a- model, ( |1.1|) . The basic 
idea is to investigate the system at large particle density when both perturbation 
theory in the coupling A and a calculation based on the S'-matrix of the massive 
physical particles apply. (For a review, see Ref. 0].) To achieve this one introduces 
a chemical potential, h, coupled to a Noether charge of the corresponding symmetry, 
and considers the ground state of the modified Hamiltonian 

H = Ho-hQ. (1.5) 

The Legendre transform of the ground state energy density or (zero temperature) 
free energy of the system must be the form: 

^ ^ 5fih) ^ f{h) - /(O) = -h^F,{h/m, Q) , (1.6) 



where m is the mass of the particles. (Eq. ( |1.6| ) follows simply from dimensional 
analysis.) The Thermodynamical Bethe Ansatz (TBA) method based on the S- 
matrix leads to a set of integral equations, which can be solved, for large chemical 
potential h, in the form of an asymptotic series expansion in h/m using a modified 
Wiener-Hopf technique |@] . 

On the other hand in an asymptotically free theory the result of a perturbative 
computation of Fq (which is of course renormalization group invariant) is an asymp- 
totic series in the running coupling, X{h/A) ex [ln(/i/A)]~^, where the A parameter 
(of dimension mass) is the usual renormalization group invariant combination of the 
ultraviolet cutoff and the bare coupling, A. 

Since JF is a physical quantity, the results obtained by TBA method and by 
perturbation theory should agree. Thus comparison of the two results provides a 
rather stringent consistency check on the S'-matrix (and also on the self-consistency 
of the hypothesis used in the course of the calculation), moreover one can also 
extract the exact m/A ratio, a rather important non-perturbative parameter of the 
theory which can then be also measured e.g. by lattice simulations. This comparison, 
following the pioneering work in Refs. [^ |l^, has since been applied to many bosonic. 



fermionic and supersymmetric integrable models. 

In Ref. a two-parameter family of integrable models has been studied by 
similar techniques. The models were mainly treated as perturbed conformal field 
theories allowing for a proof of integrability at the quantum level. The SU(2) xU(l) 
symmetric S'-matrix ( |1.4| ) is a limiting case, [pi —^ cxd, p2 = p), of the S'-matrix 
of Ref. [^. In the present paper, however, we put the emphasis on the a-model 
representation suitable for traditional weak-coupling perturbation theory. Although 
the model ( p..lD corresponds to a limiting case of the two-parameter family of Ref. , 
one cannot read off directly the asymptotic expansions for the ground state energies 
from Ref. 0. This is due to the fact that both in perturbation theory and using 
the TBA only the asymptotic series are computable. The coefficients of these series 
diverge in the pi — * oo limit, thus one cannot directly recover the coefficients of the 
limiting case from those valid for finite values of the parameters. For this reason we 
had to perform the TBA analysis directly in the p —> oo limit. 

As there are two conserved Noether charges in the model (|1.1|), Ql, Qr, there 
are three types of free energies corresponding to ground states either coupled to 
Ql or to Qr or to both of them. In fact when the ground state is coupled to Ql, 
one has to distinguish between two different cases namely p > 1 and p < 1 as 
then bound states (breathers) appear in the spectrum. We have calculated all the 
(four) different types of (zero temperature) free energies both by the TBA method 
and by perturbation theory in the modified minimal substraction scheme (MS) and 
extracted the corresponding (four) m/A ratios. Three of the kernels of the integral 
equation determining the ground state energies are simply the pi ^ oo limits of 
the corresponding ones in H] whereas in the p < 1 case with bound states we have 



computed the corresponding kernel directly using the bootstrap fusion method. We 
have found complete consistency for all four cases, in particular all of the four m/ K 
ratios agree, the result being: 



fn 3„£ p_i r(i 



- ^et-^^^, (1.7) 

where m is the mass of the SU(2) doublets. The p ^ limit (corresponding to the 
running deformation parameter g{h) —>■ —1 in (|1.1|) is expected to give the 0(3) 
NLS model (and a decoupled free field) 0. Our results are in perfect agreement 
with his expectation, as the previously calculated rrii/Ay^ = 8/e ratio in the 0(3) 
NLS model (where now mi denotes the mass of the triplet in the 0(3) model) is 
reproduced in the p — *> limit. This can be immediately seen from Eq. (|1.7|) using 
the fact that in the p — *> limit mi — ^ irpm. 

We have also proposed a lattice action for the deformed 0(4) NLS model and 
calculated the ratio of the A-parameters corresponding to the two regularization 
schemes: 

^ = y32e^. (1.8) 

Al 

The plan of the paper is as follows: in Section 2 the calculation of the free 
energies by the TBA method is outlined, in Section 3 the perturbative computation 
is presented in some detail and the m/A ratios are also given there. Section 4 is 
devoted to the perturbative calculation of the free energy using lattice regularization 
leading to the Ajj^/A^ ratio. To make the paper rather self contained we included 
four Appendices with some details of the computations. Appendix A contains some 
useful formulae to obtain the free energy by solving the pertinent integral equation. 
In Appendix B the scattering phase shift of the lowest mass bound state is calculated. 
Appendix C contains some details of the perturbative calculations, while Appendix 
D is devoted to renormalization group considerations. 

2. Calculation of the free energies from the integral equation 



As the model ( |1.1| ) admits G = SU(2)LxU(l)Fi, as its (global) symmetry group, 
one can introduce two chemical potentials coupled to the two Noether charges, 
Ql, Qr, corresponding to the two U(l) (L and R) subgroups of G. The modified 
Hamiltonian then reads 

H = HQ-hj^Qj^-hj,Qn. (2.1) 

For definiteness we chose /^L; ^r > and normalize the charges of the four particles 
as 

(1,1), (-1,1), (1,-1) (-1,-1). (2.2) 

We distinguish between the three different types of finite density ground states 
depending on hi^,h^: 



1. /iL,/iR>0, (DIAG) 

2. /iL = , /iR = /i > , (RIGHT) 

3. /iL = /i > , /iR = . (LEFT) 

In the following we shall analyse these three cases. 

2.1 The diagonal current 

We recall that our strategy will be to calculate the ground state energy density 
of the system at large particle density (and at zero temperature), that is we consider 
the system with hi^/m, h-^/m ^ 1 (m is the mass of the particles). Then it is clear 
that particles of charge (1,1) condense into the vacuum. It is less clear what other 
kind of particles (necessarily with a smaller charge/mass ratio) will appear in the 
vacuum state. We shall assume that the vacuum consist of only particles of charge 
(1,1). This seemingly radical assumption has apparently worked in all analogous 
examples studied so far, and it greatly simplifies the solution of TEA equations. 

The scattering phase of the (1,1) type particle, 5(6'), can be read off from Eq. 
( [B.ll| ) in Appendix B: 

6{e) = 6^{e) + 6^{e) . (2.3) 

When the ground state is assumed to contain only particles of charge (1,1) the 
calculation of the free energy using the TBA method reduces to the solution of 
a single integral equation ( A.l ) where the kernel K is given by the logarithmic 



derivative of the relevant ^-matrix element. 

In the present case S{6) = exp {i6{9)} where scattering phase, S{9), is given by 
Eq. (|2.3| ). We note that the effective chemical potential, h, in Eq. ( |A. 1| ) is just 



h = hi^ + /iR. From Eqs. (|2.3| , p.8| ) the Fourier transform of the relevant kernel is 

easily found to be: 

tann -rr pttIwi 



with p > 0. The result of the splitting of 1 — Ki){uj) = [G^{uj)G^{uj)]~^ where 
G+(u;) resp. G^{u!) are analytic functions in the upper resp. lower half plane can be 
written as 



2/1 iio^ 



T'(l- 



G^co) = 7^— -^-AiZ^l-_eT(i-''^(-T))e-i'^2 (2.6) 



and GJu) =G,(-u) 



For the present (fermionic type) kernel ( |2.5[ ) the constants in Eqs. ( |A.5| - |A.6| ) of 
Appendix A are given as 



a = — 



k 



2 ' 



b = Plln^- 



l-r'(l) +ln2 



(2.7) 



2p. 



GJi] 



TT 



v{?P_g-f(l-lnf) 



2 1(1 + I 



Then the expression for free energy ( |A.5| ) finally reads: 

P 






_ _ p^hl + ^ + o(hl^ 



(2. 



where t = ln{h/m) and the constant A, which enters in our final result for the m/A 
ratio is given as 



A = ^- 



Inr(l + |) + |(^-ln2)-l + 31n2 



Invr 



(2.9) 



We remark here that from the classical term in Eq. ( |2.8| ), —h'^p/n, one can imme- 
diately read off the level of an underlying ultraviolet (UV) current algebra, k, as 

k = p in|]. 



2.2 The right current 

In analogy to the assumption made in the diagonal case we now expect the 
ground state to consist of a mixture of particles of the same charge with respect 
to U(1)r, i.e. those of charge (1,1) and (-1,1) with equal densities. In Ref. [|] the 
result of the diagonalisation of the pertinent coupled TBA system has been given 
and the relevant kernel, Kii{uj) turns out to be the pi ^ oo limit in Eq. (57) of 
Ref. 0. From this one finds (somewhat surprisingly) K^{u) = Ky){u!), implying 
that Sfii{h) = Sf^ih), where we have denoted the chemical potential for the right 
current as /ir = h. So the free energies are precisely the same as functions of the 
effective chemical potentials in these two apparently rather different cases. We do 
not have an explanation of this fact, we note, however, that as long as pi < oo 

2.3 The left current 

Analogously to the RIGHT case the ground state is expected to consist of an 
equal density mixture of particles of the same charge, this time with respect to 
U(1)l, i.e. those of charge (1,1) and (1,-1). Using again the result of Ref. the 
Fourier transform of the relevant (bosonic type) kernel is found to be: 



Kl{uj) 



tanh ^^ 
2 cosh ^ 



(2.10) 



where x = p — 1 and x is assumed to be positive. The functions appearing in the 
decomposition of 1 — Kl{uj) are now given as 

r (fA- ^^ r(i - f ) i;^fi_in(-H^)), .^-.-.s 

and G_(co') = (?+(— cj). For the kernel ( |2.1CI| ) the constants in Eq. ( |A.8| ) are not 
difficult to calculate and we find: 

a = -f , b= — [xlnx + x(ln2-l -r(l)) + 21n2] , 

k = ^, G+« = ^;^e-f(^-"^f). (2.12) 

So the asymptotic series of the free energy density can be simply written as 

^hW = -^ (^t + (1 - |) int + A + . . .) , (2.13) 

where the constant A, which enters to our final result for the ratio m/A is found to 
be 

A = -ln2-ln7r-l + lnrfi^) +-(l-ln2). (2.14) 



At this point we would like to remark that although the p — > oo limit in the 
S"- matrix is smooth (c.f. Eq. (p.8| )) this is not the case in the final results for the 
free energy, i.e. in Eqs. (|27^), ( |2.13| ). This is not so surprising since Eqs. (p7^), ( |2.13| ) 
have been obtained as asymptotic expansions valid for h/m — * oo for fixed p, and 
the lack of a smooth p —>■ oo limit of the coefficients shows that the h/m —>■ cxd and 
the p — > cxD limits cannot be simply interchanged in these asymptotic expansions. 

Next it might be worth pointing out the following (somewhat puzzling) fact. In 
Eq. ( |2.8|) the ratio of the coefficients of the Int/t^ term, — p^/4, and the — 1/t term, 
p/2, is Td = —p/2, while in Eq. ( |2.13| ) the ratio of the coefficients of the Int term, 
1 —p/2, and that of the t term is Tl = 1 —p/2. This appears to be a serious problem 
since on general grounds in a theory with a coupling. A, the perturbative result for 
a fermionic type free energy is expected to be of the form 

ShAh) = -^ (l - /i^ + OCX')] , (2.15) 

while in the same theory for a bosonic case it is expected to be given by 

^hoAh) = -^ (" J- + const. + 0{X)] , (2.16) 

where /o, /i, 6o are constants and X{h) is the running coupling. Since from the 
renormalization group 

^^A.,4 + |.„.,4 + 0(J^). (2,17) 
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one would expect rper = '"Bos = /9i//5o- As in the present case r^ 7^ Tl, it seems 
impossible to match Eq. ( |2.^ ) and Eq. ( |2.13| ) simultaneously with the corresponding 
perturbative expansion. Fortunately this turns out to be only an apparent paradox 
which is resolved, however, in a not entirely obvious way (see Section 3). 

Finally we note that the leading term of Sfi^{h)/h'^ diverges as \n{h/m) which 
can be interpreted that either the underlying UV current algebra does not exist or 
that if it could be defined its level is divergent. 

2.4 The left current for p < 1 

In the previous subsection it has been assumed that the parameter, x = p—1 > 0. 
Although in the S-matrix any p > is allowed in Eq. ( |2.1CI| ) the range of p is. 



however, restricted to p > 1. This restriction follows from the fact that in Eq. ( 2.10 ) 
lim^^oo -^l(^) = 1 for < p < 1, that is for p < 1 Ki^{uj) ceases to be the Fourier 
transform of a function. Then clearly one cannot just simply analytically continue 
in p the result for the free energy, (|2.13|) , for p < 1. 

In order to clarify what happens when p becomes smaller than one, we recall 
that then in the Sine-Gordon (SG) theory bound states (breathers) appear in the 
spectrum (for p < 1 one enters simply the attractive regime). The SG mass spectrum 
is given by the well known formula: 

TTpr 1 

m^ = 2msin-^, r = l,2...<-, (2.18) 

2 p 

where m denotes the mass of the SG kinks. As the S matrix of the present model, Eq. 
( |1.4D , contains the SG 5* matrix as a factor the spectrum then changes accordingly, 
i.e. breathers also appear (with apropriate multiplicity) in the deformed NLS model. 
(We remark that m is also the mass of the SU(2) doublets, of course). 

The appearence of these new particles implies that as soon as p becomes smaller 
than 1, the vacuum changes as then it becomes energetically favourable for the r = 1 
charge (2,0) breather to condense. This charge (2,0) breather can be interpreted as 
the bound state of the (1,1) and (1,-1) particles (kinks). In fact as this particle 
has the highest charge/mass ratio, we can again assume (as in Subsection 2.1) that 
the true vacuum consists only of the condensate of the r = 1 charge (2,0) breather. 
Then one needs the scattering phase, 5^^''^\9), of the (2,0) breathers in order to 
compute Sfi^{h) for the p < 1 case. 

The calculation of 5^"^'^^ (9) amounts to a simple application of the well known 
bootstrap- fusion method |T3[. In Appendix B we summarized this method and 
calculated the breather-breather phase shift, which is given by Eq. ( p.l6|) . Using 
Eq. (|B.18| ) the Fourier transform of the pertinent kernel is found to be 

1 -irB(cj) = 2tanh^^cosh —e —, (2.19) 



where |a;| = |p— 1| = 1 — p. The functions appearing in the decomposition of 

1 — Kb{uj) are now given as 






G+(^) = -^=^ ' 'n i.. ' e-'^^^-'^i-^)) , (2.20) 



and again G^{uj) = G'+(— cj). 

For the kernel ( p.l9 ) the constants in (|A.8|) are found to be: 



a= ^, 6 = -[|a;|ln|a;| + |x|(ln2-l-r(l))-21n2] , 

k = ^, G40 = 2r(^-^je^(^-^°y). (2.21) 

So the final result for the asymptotic expansion of the free energy density for the 
r = 1 charge (2,0) breather condensate can be written as 



5/b(M = -^ (r+(l-|)lnr + AB + ... 



f2.22) 



where r = Inh-B/mi and the constant Ab is given as 

AB = 31n2-^-ln7r + lnr(^0+|(l-ln2)+ln(^sin^) . (2.23) 

It is now easy to verify that expressing the breather mass, mi, through the kink 
mass, m, using the SG spectrum formula Eq. (|2.18| ) and taking into account that 
the charge of the breather is twice that of the kink, i.e. Hb = 2h, the free energy 
density of the breather condensate, SfB{hB,mi), becomes identical to that of the 
charge (1, ±1) kink mixture, 6fi^{h,m) analytically continued for negative values of 
X. This is very reassuring as the result for Sf-L{h, A) from perturbation theory turns 
out to be continuous for the full range < p < oo, which is certainly to be expected. 

2.5 The 0(3) model as the p — » limit. 

In this last subsection we discuss the p — > limit and recapitulate some argu- 
ments advanced a long time ago by Wiegmann [Q that the SU(2)lxU(1)r invariant 
theory ( |1 . 1|) in this limit becomes equivalent to the 0(3) a-model and a decoupled 
free massive particle. First of all from Eq. (|2.18| ) one sees that m/rrii — > cxd as p ^ 
i.e. the kinks 'disappear' from the spectrum. Furthermore the higher breathers with 
r > 1 also 'disappear' from the spectrum, but for a different reason. The bind- 
ing energy of the r > 1 breathers (which can be considered as bound states of the 
r = 1 ones) tends to zero giving rise to zero energy bound states. Therefore the 
limiting spectrum consists of 4 massive particles, transforming as (3l + 1l, 1r) under 
SU(2)lxU(1)r. As already claimed in Ref. and confirmed in Ref. [^ in the p — > 
limit the S-matrix of these remaining r = 1 kinks becomes effectively 5*0(3) ® 1 i-e. 
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the tensor product of the (minimal) 5* matrix of the 0(3) model proposed by the 
Zamolodchikov brothers ||T^ with the S-matrix of a free boson, consistently with 
this somewhat heuristic argument. As this is a highly non-trivial limit it is clearly 
desirable to perform as many quantitative checks on it as possible. In the following 
we add yet another piece of evidence corroborating the arguments of Wiegmann, 
by showing that the p ^ limit of the rrii/Ajj^ ratio agrees precisely with the one 
calculated directly in the 0(3) model 0. 

3. Perturbative Calculations 

In this section we calculate the ground state energy of the system in the pres- 
ence of external fields in perturbation theory. As in the previous section, we shall 
consider fields coupled to the Noether charges corresponding to the U(1)lxU(1)r 
transformation 

SG = leLcr^G - ienCa^ . (3.1) 

To calculate the ground state energy one has to gauge the symmetries in the Eu- 
clidean field theory formulation and consider the case of constant, imaginary gauge 
fields 0. The corresponding covariant derivative D^ is given by 

D2G = diG + hLa^G - hnGa^ , D^G = d^G . (3.2) 

Putting ( t3.2|) into ( |1.1D we write the gauged Lagrangian as: 

£ = £0 + >Ci + £2 , (3.3) 

where £1 and £2 denote the terms linear and quadratic in the external fields, re- 
spectively. 

To actually calculate the free energy in the presence of the external field, we take 
the master formula 

Hh) = Sfih) = fih) - /(O) , (3.4) 

where 

^-Jd"xfih) ^ fvGexp[- /'t/"x(£o + £i + £2)}, (3.5) 

and expand it in powers of the (bare) coupling Aq. From now on the bare lagrangian 
parameters (couplings) will be denoted as Aq and go. We have also written the 
volume element as d'^x, where n = 2 — e to indicate that we are going to use 
dimensional regularization in the perturbative calculations. 

In perturbation theory the action has to be expanded around a (stable) solution 
of the classical equations of motion. Thus our next task is to find elements of 
SU(2), Go, that correspond to (local) minima of the quadratic part of the gauged 
Lagrangian, 

2hl 2^ 24(1 + (7o) , 4hLhR{l+go)z 

1^2 = r— (1+^0^) T^ + 7 , (3.6) 

Aq Aq Aq 



where 

z = ^Tr {a^G- V-^G} , (3.7) 

which plays the role of the potential energy. In this paper we shall consider the 
following two solutions. The first one, called 'bosonic', (BOS) is given as: 

Go = -^(l + za^), hn = 0. (3.8) 

The other one, which shall be referred to as 'fermionic' (FER) is written as 

Go = i(^^ , (3.9) 

where we also require that h^ and hji satisfy the inequality 

A = A[hlgo + hLhnil + go)] > , (3.10) 

ensuring the non-negativity of the mass term. 

The terms 'bosonic' (resp. 'fermionic') are used here by anticipating that they 
actually correspond to the bosonic (resp. fermionic) type free energies discussed in 
the previous section. The BOS-type solution will be used to calculate the free energy 
corresponding to the pure U(1)l charge (/ir = 0, LEFT) case, while the FER-type 
solution is appropriate to calculate the free energy for both the pure U(1)r charge 
{Hl = 0, RIGHT) and the 'diagonal' (DIAG) case, where both hi and Hr are 
different from zero. 

We now start to analyse the BOS case first. The leading, 0{Xq^), term in 
perturbation theory is given by the potential energy, C2, at its minimum, which is 
immediately seen to be {zq = 0, hji = 0) 

^(-i) = /:(-i) = _^. (3.11) 

Then we turn to the calculation of the next, 1-loop, term in perturbation theory. 
The calculation is elementary but for completeness we give some details in Appendix 
C The general 1-loop contribution is given by Eq. ( |C15| ), which is a rather compli- 



cated function of three parameters. It simplifies somewhat in the two cases we are 
interested in. In the present, BOS, case after using ( |U.12|) , we obtain 



WO) ^ ^ / JH!. {^+90)P2-90P /o 10X 



4/ii f _d^ {'^ + go) pI - gop"^ 

n J (27r)" (p2)2 _ 4g^hlp^ + 4(1 + g^)hlpl 

Because of the explicit dependence of the integrand on the second component of the 
momentum we have to modify the usual rules of dimensional regularization by the 
following. In all Feynman integrals we perform the one- dimensional integral over 
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P2 first, followed by the integration over the remaining momentum components, 
continued to n — 1 dimensions, i.e. in the present paper 



(27r)" ■ J (27r)"-i J (27r) ' 



-CXD 



It is possible to calculate the integral ( p.l2| ) exactly |]T6[. For our purposes, however 



it is sufficient to calculate its divergent part (the part that becomes singular in the 
e ^ limit) and to show that JF(°) must be of the form 

J'^'^ = hlU^ + Wigo)y (3.14) 

where W{gQ) is a (complicated) finite function of its argument. Fortunately all we 
shall need is the value of W{go) at the special point go = —1. For this very special 
go ( ^.121 ) simplifies enormously: 



^(°) - ^^i / ^> ^ (3 15) 

9o=-i n J (27r)" p2 + 4/^2 • l'^-^^^ 

Calculating (|3.15|) and comparing the result to (|3.14|) yields 

Vr(-l) = ;^[l + r(l) + ln7r]. (3.16) 

ZTT 

Having obtained the bare 1-loop result, we now turn to renormalization and renor- 
malization group (RG) improvements. The methods are well-known but we sum- 
marized the details of the calculation in Appendix D to make this paper relatively 
self-contained. 

The RG improved perturbative result gives the asymptotic expansion of the 
free energy for large values of the external fields. Writing the sum of the classical 
contribution (|3.11|) and the one-loop correction (|3.14| ) and using definitions ( p. 2]) , 



( p.3| ) we obtain the renormalized form 

^(/,) = -^ _ (1^^ In ^ + h'w{g) + 0{X) , (3.17) 



where h = hi. Expressing the result in terms of the single effective coupling, Ach, 
defined by Eq. dUJoD we get 

^,„)_2„.{_L + _L(,„2-i) + ...j. (3.18) 

At this point we remark that we have identified the RG invariant parameter defined 
in terms of the running couplings g{t) and A(t) by Eq. ( p.21|) with the 'deformation 
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parameter' p appearing in the S matrix ( |1.4| ). Finally expanding ( |3.18| ) with the 
help of ( |D.32| ) the free energy in the BOS case can be expressed as 



^(ft) = _^j,+ (l_|),„,+ (,n2-i) + ...j. 



(3.19) 



Next we turn to the computation in the FER case. The leading (classical) term 
of the free energy is then 

^M) = ^(-1) = _^(i+M ^h, + h^f . (3.20) 

Note that JF*-^^-* for the diagonal charge depends only on the sum hi + Hr. This is 
consistent with the results of the TBA analysis found in Section 2. for the DIAG 
case. In the present (FER) case ( p.l5|) becomes 



^ ~ nJ (27r)" '^%2 + A - 2tcup2 ' ^ ^ 

where 

u = hiil - go) - hnil + go) . (3.22) 

Evaluating Eq. ( |3.21| ) using the modified dimensional scheme (c.f. ( |3.13|) ) yields: 

^(°) = (1 + go)"" {hi + M"|^ + i [^ + ^'(^) + ^^H I • (3.23) 

It is gratifying to find that the one loop term (|3.23|) also depends only on the sum 
hi + hR, just like the classical one ( |3.20| ). This is not an accident. It is in fact 
absolutely essential that J-'^hi, h^) be a function of hi + h^ in order to be able to 
match the perturbative result with the one found by the TBA method for the DIAG 
case in Subsection 2.1. As shown in Appendix C, in the FER case the free energy 
does depend only on the sum h^ + h^ to all orders in perturbation theory. Moreover 
it is easy to see that the free energy in the FER case must be of the form 

^Y^ f: FLigo) (XoM-^)"^] , (3.25) 



-2/.^ i±^ 



where 

M={l+go){hL + hn), (3.26) 

and the summation is over the number of loops, L. (|3.24D is nothing but dimensional 



analysis based on the fact that the only dimensionful parameter available is M, which 
serves as an infrared cutoff in our perturbative computations. The second form 
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( |3.25[ ) illustrates the fact that all higher order terms are proportional to (1 + goY, 
which fact plays an important role in our calculations. 

Adding now the classical and the one-loop terms, ( |3.2CI| ) and ( |3.23| ) and expressing 
the result in terms of renormalized quantities we obtain: 




2U2' 



In 



[i + gn 



1 - 27 + 0(A) 



where h = h^ + hji and the constant 7 is defined as 



7 



ir(l) + lnV47r. 



(3.27) 



(3.28) 



As we need the expansion of J-'{h) up to 0(Afg) in the effective coupling in this case, 
this would necessitate a three-loop computation at first sight. In fact we obtain 
all the 0(Afg) terms from the one loop result ( p. 27 ) alone! This 'mini miracle' is 
due to the {1 + g)"^ factor in front of the one- loop term together with the fact that 
(1 + g)^ = C(Af(j), implying that the one-loop term is already C(Afjj.). Moreover as 
it follows from Eq. ( |3.25| ) the contribution of the higher order terms will be at least 
O(Afg) hence they can be safely ignored. Putting then everything together we find 



Hh) 



ph"^ 



n 



P A + 



327r2 °ff ^ 87r2 



A. 



In p + In ( ^— 



0{Xl)\. (3.29) 



Note the non-analytic contribution, InA^ff, in the last term of (|3.29|) . This comes 

ultimately from the (^o-dependence in ( p.26| ) of M, which plays the role of the infrared 

cutoff in our calculations. The presence of this non-analytic term in X^g explains 

why the coefficient of the Ins term differs for the BOS and the FER cases. Indeed 

using Eq. ( |D.32| ) we find that the large s expansion in the FER case is finally given 

as 

,2 



Hh) 



ph 



TT 



p_ 

2s 



— ^ In s + — ^ 
4s2 2s2 



Inp + 



P 



+ ■ 



(3.30) 



Having calculated the asymptotic form of the free energy both in perturbation 
theory and with the TBA method, by comparing them (with t = s — ln(r?7,/Af^)) 
we obtain the relation between the mass of the doublet particles, m, and Aj^. In 
both the BOS resp. FER case the comparison of Eqs. (|2.13|) and (|3.19| ) resp. 
and ( p.30| ) leads to the result ( |OD as already announced in the Introduction. 



We emphasize that it is already a very nontrivial check on the overall consistency 
of our assumptions that all the expansion coefficients, for both the BOS and the 
FER cases, agree. A further, very stringent consistency check is that the m/A ratio 
obtained in the BOS case is exactly the same as the one obtained in the FER case. 
Finally by comparing the free energy of the breather condensate, computed from 
the TBA method, Eq. (|2.22| ), with the perturbative one, Eq. ( p.l9|) , one finds 



mi 



A 



E E_1 
"2 62 



Ffl 



MS 



p 2sinf 
2 np 



(3.31) 
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By recalling that the mass of the lowest lying (r = 1) breather is mi = 2m sin ^ one 
immediately obtains from Eq. ( p.31|) the r/i/A ratio in Eq. (|1.7| ). This is clearly a fur- 



ther nontrivial check on our result and on the mutual consistency of the hypothesises 
made in the course of the calculation. 

The ratio ( |1.7D diverges in the p — > oo limit. This illustrates the fact, discussed 
in Subsection 2.3, that it is not possible to recover the asymptotic expansion of 
the ground state energy of the 0(4) model from the generic case in this limit. 
Furthermore m/Ajj^ also diverges as p — > consistently with the Wiegmann scenario 
as discussed in Subsection 2.5, while 

,. m^i 8 
lim 



which exactly reproduces the m/A ratio of the 0(3) model first computed in Ref. 
0. In our view this beautiful agreement certainly yields some more quantitative 
support as to the correctness of the 'conventional wisdom'. 

4. Lattice regularization 

In this section we propose a lattice action for the deformed principal model and 
calculate the free energy in this lattice version of the model. This is then used to 
calculate the ratio of the A-parameters of the dimensionally regulated theory and 
the lattice version. 

First we rewrite the Lagrangian ( |1 . 1| ) in the form: 

Co = {l + g)C^^'^-gC^^'\ (4.1) 

exhibiting that it interpolates between the 0(3) and 0(4) nonlinear cr-models. In 
fact 

CO(^) = ^tr{d,G-%G} , (4.2) 

is the familiar action of the SU(2) principal chiral (or 0(4)) nonlinear a- model and 

£^(3)= l5^S'^5^S% (4.3) 

is the Lagrangian of an 0(3) symmetric nonlinear a-model with a composite 0(3) 
field 

S'^ = itr{G"V"G(T^}. (4.4) 

Motivated by this form of the Lagrangian we take as our lattice action the 
following linear combination. 

S,.u-...= {l + ~9o)S'"~^^-goS'''-'K (4.5) 
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Here 

50(4) = i ^ tr {2 - G-\x)G{x + fi)- G-\x + ix)G{x)} , (4.6) 

is the standard 0(4) lattice action in terms of the SU(2) principal model variable 
G{x) associated to the lattice site x. (As usual, the lattice spacing a is taken to be 
unity and fi denotes the unit vector in the fi direction. Moreover, we have denoted 
the bare lattice couplings by Aq and ^0 to distinguish them from their counterparts 
in the dimensional scheme.) Similarly 

gois) ^ i^ |i_2»(x)S"(x + /i)} (4.7) 



A 



x,^ 



i 5] tr {G{x)a^G-\x)G{x + fi)a^G-\x + /i) - l} , 
2Ao x,fi 



is the standard 0(3) lattice action in terms of the composite field (|4.4| ). 



For simplicity we shall consider here the free energy associated to the U{l)i^ 
charge only. This corresponds to the BOS case of Section 3 (i.e. hi = h, hji = 0). 
An apropriate gauging of the Lagrangian ( |4.5| ) and taking imaginary gauge fields 
leads to: 

5,a.,ed = ^^Y.^^{^-G-\x)e~'^-'^'G{x + fi)-G-\x + fi)e'^-'''G{x)} 
Ao x,fj. 



+ -^ ^ tr {G{x)a^G'\x)e-''^'''G{x + fi)a''G-\x + /i)e'^-'^' - l} 
2Ao x,ii 



with h^ = h6f,2- 

The perturbative calculation of the free energy on the lattice is now completely 
analogous to the one in the dimensional scheme. The leading contribution is given 
again by the minimum of the gauged action for constant fields. (|4.8| ) also takes its 
minimum at the group element (|3.8|) and is given by 



4-^) = -lii±M Uoshha-l) + J^ (cosh2/.a- 1) , (4.9) 

where the a-dependence has been reconstructed. In the continuum limit, a — ;> 0, 
(IJ) reproduces ( pl| ). 



The 1-loop contribution is half the logarithmic determinant of the operator cor- 
responding to the quadratic term in the expansion of the lattice action around the 
group element ( |3.8| ). We have not tried to calculate the 1-loop determinant for 
generic ^o; because for our purposes it is sufficient to show that it must be of the 
form 

^f = -^'{^^ In/ia + W^Cgo) + 0{a')\ . (4.10) 
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PVl(^o) similarly to its dimensional counterpart, is a complicated function of its 
argument but all we are going to need is its value at the 0(3) point qq = —1, where 
it is relatively easy to evaluate. We find 

W,{-l) = -^-^\n2-l. (4.11) 

ZTT ZTT O 

We now proceed, as before, from the bare results (|4.9|) and ( [4.10| ) via the renor- 
malized ones to the RG improved 1-loop formula 

-i-2--fU.... (4.12, 




where we have used the results of Appendix D. Comparing ( ^.191 ) and ( [4.12] ) gives 
the ratio of the A-parameters of the two regularization schemes: 



Ar 



57r 



32 e". (4.13) 



Note that the ratio ( 4.13|) is independent of the RG- invariant parameter p, but 



is nevertheless different from the corresponding ratio of the standard 0(3) model, 
which is given by a formula similar to ( [4.13|) , but with an exponent that is 7r/2 
instead of Svr/S. This means that our lattice regularization is different from the 
standard one in the 0(3) limit. A possible explanation for this might have been 
that the difference is due to the compositeness of our 0(3) field (^4.41). But by a 



simple computation we find the standard A^^ for qq = —1 ii the first term in (|4.5|) is 
simply omitted. To clarify this apparent paradox we note that the difference is due 
to the following term in the quadratic piece of the gauged lattice action. 

S^ = '^^'^ + 9o) ^cosh/i^$(x)[<l>(x)-<l>(x + /i)l. (4.14) 



Ao 



X,IJ. 



( [4.14| ) is, of course, absent if ^o = —1- This is the case in the standard lattice version 
of the 0(3) model. If, however, ^o ¥" "1; then the coupling dependence can be scaled 
out of ( |4.14] ). Close to the continuum limit, after rescaling, ( |4.14| ) corresponds to 
the continuum Lagrangian 

C^ = l d^(t> d^ct> + ^ h'^a' {d2(t>f + 0(/iV) , (4.15) 

from which one would naively conclude that this field is decoupled from the external 
field h in the continuum limit. This is, however, not correct since the second term 
in ( |4.15| ) does in fact contribute to the continuum free energy, because of quadratic 



divergences in lattice perturbation theory. This contribution is /i^/8 leading to the 
modified exponent in ( [4.13|) . We expect to find the same modified A^/A^jg ratio by 
comparing the 4-point function computed in the two regularization schemes. 
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For completeness we also computed the A^jg/A^ ratio in the 0(4) symmetric 
case g = 0. In this special case the function defined in (|3.14|) is relatively easy to 
calculate. In the dimensional regularization scheme 

iy(0) = l±^, (4.16) 

whereas in the analogous calculation for lattice regularization we find 



1 



TT 



An 
Using ( [4.16| ) and ( |4.17| ) we obtain the Aj^/A^, ratio 

Ar 



VrL(0) = — 1 + In32 + -J. (4.17) 



32 64, (4.18) 



in agreement with the well-known result. We note that what is new in our calculation 
is that we considered the external field coupled to Ql here instead of the diagonal 
charge considered previously. The A^j^/Ai ratio is the same of course, as it should. 

We think that it would be extremely interesting to perform Monte Carlo simu- 



lations using the lattice action (O). This would provide us with an independent 
way of determining the ratio of the mass of the physical particles and the lattice A 
parameter and also with an interesting example of a lattice field theory model with 
more than one relevant couplings. 

5. Summary and Conclusions 

In this paper we have computed the ground state energies of the deformed prin- 
cipal model in the presence of external fields for all three possible orientations of the 
external field. We have found a complete consistency between perturbation theory 
and the results from the TBA method based on the proposed S-matrix of the model. 
In all three cases the leading and the first few subleading coefficients of the large 
field expansion for the ground state energy are identical once we calculate the m/A 
ratio using one of the expansions. This beautiful agreement is absolutely nontrivial 
since it can only be achieved after identifying an unusual non-analytic term in the 
perturbative expansion for the fermionic case. Such non-analytic terms can only 
occur in perturbation theory in the presence of several couplings. To our knowledge 
our example is the first one where the results of the TBA program is compared to the 
results of asymptotically free perturbation theory with several coupling constants. 

Actually the presence of the running deformation parameter, g, made the compu- 
tations much easier as compared to a purely fermionic model. In a fermionic model, 
generically a three-loop calculation is necessary to obtain the m/A ratio, whereas in 
the deformed model here, it has been sufficient to calculate only to one-loop order! 
This nice simplification is due to the fact that all higher terms are proportional to 
{g + 1)^, and since for large fields g -^ —1, they are effectively supressed. 
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We have also outlined a proof to all orders of perturbation theory that the 
free energies for the two apparently different fermionic cases (DIAG and RIGHT) 
actually have the same form as functions of the effective chemical potentials. We 
have no explanation for this fact but it is yet another very nontrivial consistency 
check since the TEA integral equations are identical for the two cases. 

It is also quite remarkable that the agreement between perturbation theory and 
the TEA method is valid for all values of the RG invariant parameter, p = 2tt{1 + 
g)/X, including the attractive regime < p < 1, where the bootstrap S-matrix 
indicates the existence of a complex bound state structure in the spectrum of the 
model. We have found that the structure of the ground state is changing with a 
decreasing p in the sense that it is always the lowest lying breather bound state 
that condenses to the vacuum. The final formula for the ground state energy is, 
however, analytic for the whole range of < p < oo and agrees with the results of 
PT. Finally our results are also in perfect agreement with the Wiegmann scenario 
for the p — i> limit since linip^o ^i/-^ = 8/e, {rrii being the mass of the lowest lying 
breather), thus we reproduce the m/A ratio computed previously in the 0(3) model. 
What we found here lends also additional support to the empirical rule that only the 
particle with the highest charge/mass ratio condenses in the ground state. This rule 
was found to be correct in all cases considered so far and proved to be correct 
also for the deformed model in the attractive regime. In view of the complexity of 
the spectrum of physical exitations in the attractive regime it would be extremely 
interesting to find an independent explanation of this phenomenon. 

We have also proposed a simple lattice version of the SU(2)xU(l) symmetric 
model and computed the Aj^/A^ ratio where A^ denotes the lattice lambda pa- 
rameter. 

An interesting question that requires further study is the application of the 
TEA program to calculate thermodynamical quantities for finite temperature and 
comparing them to the results of perturbative calculations. 

Acknowledgements J. E. gratefully acknowledges a CNRS grant. This investigation 
was supported in part by the Hungarian National Science Fund (OTKA) under 
T019917 and T 030099. 

Appendix A. Solution of the TBA integral equation 

In the Appendix of Ref. |TD[ the solution of the pertinent integral equation de- 



termining the free energy has been obtained for a generic kernel by applying the 
generalized Wiener-Hopf technique. The integral equation is given as follows: 

e(e)- f de'K(e-e')e(e') = h-mcoshe, (a.i) 

J-B 

together with the boundary condition e(±i?) = 0. In terms of the solution of 
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Eq. ([A.l|) the free energy is 



7D C^ 

sf{h) = -— de cosh{e)e{e) . (a.2) 



We do not go into the details of the solution technique of Eq. ( |A.1| ) as they can be 
found in the literature (e.g. Ref. [|lOl)- Here we just present the results in a somewhat 
simplified form which could prove convenient for future applications. The crucial 
step in solving Eq. ( |A.1|) by the generalized Wiener-Hopf method is to split the 
Fourier transform of the kernel in Eq. ( [A.l| ), K{uj), as 

l-Kiiu) = [G4iu)G4iu)]-\ (A.3) 

where G+(a;) resp. G-(a;) is analytic in the upper resp. lower half plane. When 
the kernel is of the 'fermionic' type, i.e. K{uj) ^ 1, one can assume the following 
expansion of G+(z^) around ^ = for ^ > 0: 

G4tO = ~kexp{-d^\n^}{l + k + Oie)) ■ (A.4) 

Then in the limit of high density h/m ^ 1, the asymptotic series for 6f{h) can be 
obtained by a simple iterative technique. The final result, Eq. (A. 34) of Ref. |TD[, 
can be conveniently recast in the following form: 

where t = ln(/i/'m) and the constant A is given as 

A = d(l-b-\n-^] +d^(- + r{l)-\n2) . (A.6) 



For 'bosonic' type kernels, K{0) = 1 so the splitting functions G^{uj), G^{lj) 
are singular at tu = 0. The solution of the integral equation ( |A.1| ) is then more 
complicated than in the previous (FER-type) case as the corresponding Neumann 
series is not uniformly convergent and the solution cannot be obtained by simple 
iteration. The solution of Eq. ( |A.1D for a generic BOS-type kernel in the limit 



t = \n{h/m) -^ oo has been given in Ref. [12| up to order of Int/t. According to 
the results of Ref. |T2| for G+{i^) admitting the following expansion around ^ = 

(e>o) 

G+(zO = ^exp{-aeinO(l-&e + 0(a) , (A.7) 

the asymptotic series for t ^ 1 of the Legendre transform of ground state energy 
density can be expressed in the form 

5f{h) = -^e [t+{a + ]^)\nt + A + ..^ , (A.8) 
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where the constant A is given by 



A = lnp^-— -l + a(ln8-l-r'(l)). (A.9) 



Appendix B. Bootstrap fusion 

In this appendix we calculate the phase shifts corresponding to the scattering 
of two bound state particles of charge (2,0) and of lowest mass mi = 2m sin ^ 
in terms of the scattering phases of the "elementary" particles of charge (1,1) and 
(1,-1). The procedure is called bootstrap fusion and being an important part of 
the whole bootstrap program is discussed at length in the literature |l^, |18|, |20 |. 



Therefore we just apply the general framework to the special case of our interest, 
but even this could prove to be useful for other models. 

We shall assume that the S-matrix describing the scattering of the "elementary" 
particles of mass m can be decomposed as 

sffio) = Y.^'e'''Se{ewi,, (B.l) 

e 

where the orthogonal projectors satisfy 

The scattering phases may have poles in the physical sheet < ImO < tt of the form 

Se{e) ^ -^ . (B.3) 



If the pole in ( B.3 ) is associated^ with a bound state of mass 



mi, = 2171 COS — , (B.4) 

then one can use the bootstrap fusion procedure to calculate the S-matrix of the 
bound states from the S-matrix of their constituents. We assume further that the 
bound state, b, occurs only in one of the chanels (which we shall label e = b) and 
that the scattering of b is diagonal with all other (elementary or composite) particles 
from the spectrum of the model. Then the scattering of b and X (which can be one 
of the constituents or some composite particle) is given by 

a,f3,X,a,(i 



^ Bound states are not the only possibility to explain the presence of these, or higher order, 
poles. See @ |o|. 
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In the case of the Sine-Gordon model the particles are indexed by their U(l) 
charges (which can be 1 or -1) and the four scattering chanels are labeled (0, 0, ±). 



The projectors are |19 



vf = ^%=s--y.^, (B.6) 

The soliton-soliton and antisoliton-antisoliton scattering phases are 

Sjf\e) = S^''\9) = -6'^"^'^ (B.7) 

where 

°° ' ( n] sinh f ^"(P"^) ^ 

SM= J du^'-^K.iu), k,{u) = ^ ^ ^\ ' ^ j_^ . (B.8) 



2 cosh (f^) sinh (^) 



The other two eigenvalues are 



sinh 



S'fi0) = -—^j^s;f\e), (B.9) 



and 



cosh 



. 2p J 
' 9+in ' 



S^I\9) = TTii^^^^W- (B.IO) 

cosh (^^) 

The lowest breather has mass mb = 2m sin ^ = rrii corresponding to Mf, = 7r(l —p) = 
Ml and it is present in the chanel e = — only. 

Now we return to the deformed principal model. The S-matrix of the model, 
Eq. (|1.4| ), is a tensor product of two copies of the SG S-matrix discussed above with 
p = oo in the first factor. Using the {Ql, Qr) labeling of particles, ( [L.4| ) actually 
means 

^{aL,aRWL,l3R)y^) - ^CLh y^J^anPR [^ ) ■ K^-^^) 

Using a tensor product labeling E = (e^; cr) also for the scattering chanels, the 
projectors can also be written in the product form 

^{aL,aR)if3L,f3R) = ^I'^Lh "^TrPr ) (B.12) 

and similarly for $£;. Finally 

^(e.;e,)W = 5lrH^)^SW- (B-13) 

We are interested in the bound state B = (0;—), which is composed of the 
elementary particles (1,1) and (1,-1) and has lowest mass mi in this charge (2,0) 
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chanel. Applying the bootstrap fusion formula ( |B.5| ) to calculate the scattering of 
B and one of its constituents (1,7_r) gives 

SB(i,j^){e) = e^'5^(^+^)e''5-(^-^)5{^)(^) , (B.14) 

where 



,(p)^fl^ _ sinh g + z sin ^^ 
sinh 9 — i sin ^' ^^^ 



SriO) = . _ . . .,j^ , (B.15) 

2 

is the analogous soliton-breather scattering in the SG model. We can interpret 
( [B.14| ) as the product of this SG scattering coming from the second factor in ( |1.4D 
multiplied by the diagonal soliton-soliton scattering coming from the first one. 



Having obtained the bound state-particle scattering, we now apply ( |B.5| ) once 
more to find the B — B scattering phase, 5*^^'°-* {9) . Using (p.l4| ) in (p.5| ) we find 

Sbb{9) = e^^'^'^'W = e^^^-We^^-^^+^i^e^-^-^'-'^i^^SV) , (B.16) 

where 

^(p),^, sinh 6^ + i sin pvr , , , , /-n^^^N 

SiiiO) = —r^ r— ^, ui = n{l-p)=7rx, B.17 

smh 9 — 1 sm pir 

is the breather-breather scattering phase in the SG model. The interpretation of 

( [B.16|) is completely analogous to that of ( |B.14|) . Finally for later convenience we 

also give here the Fourier representation of S{i{9): 

S^\9) = -exp ^ -^ [ , J^-u9 cosHl-p)nu; 

OJ cosh j^TXUJ 

Appendix C. Perturbative calculations 

We chose the following parametrization for the group valued field G: 

^ = ^°7rT^(^ 1 )('o e^O- ^^-^^ 




Here ^ is a complex scalar field, $ is a real scalar field and Gq G SU{2) is an 
appropriate constant group element given by (|3.8|) or (|37 



Before we can expand the Lagrangian in powers of the coupling constant Aq we 
have to rescale our fields as follows. 



In terms of the rescaled fields we have 



^° = (1 + 4^12)2 + 2^^^^^^^^ - 2 V ^°(^ + 9,)d,M, - ' g"^^ " AA , (C.3) 
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where 

^'~ i + ^w ■ ^ ^^ 

Due to the fact that the two cases BOS and FER correspond to expansions 
around different points, Ci and £2 are different in terms of the rescaled fields. For 
the BOS case we find 



A = V2(l + .0)/.. 7^^520 + U^JI^ (52^ -92^ 



1 + fl^P V 2 1 + fl^l^ 

Aq i^ + r 
2 i + ^w 



+ 9ohL\l- ^ . Aou/,12 -^2 



(C.5) 



and 



2hl ,2 ^ + V'% .p., 

Ao ri + ^|^i2^ 



whereas in the FER case 



1^1 



2 



Ci = 2ihL\/Xo{^ + ^0) ^20 , , Api./.i9 + 



1 + f |^|2 



/ii(l - go) - hR{l + go) + , , ^| ,, 



1 + fl 



A, 
(C.7) 



and 






The leading, O{\o ), terms of the free energy are ( p.ll| ) and ( p.20| ). The next 



terms are of the order 0{1) and are given by half the logarithmic determinant of 
the operator corresponding to the 'free' Lagrangian, that is the piece quadratic in 
the fields ip and (j). This 'free' Lagrangian can be written in the BOS case as 

£(0) = d^r d,^ + ^-9^0 ^^.ct> - gohli^ + r)' + ip{l + go)hL{^ + V) <920 , (C.9) 
and for the FER case as: 

£(0) = d^r ^^^p + ^9^0 9^0 + A|^p + uj {^d2r - rd2^) . (CIO) 

In both cases, only two fields are really coupled to the external fields. These are 
the real part of ip and in the BOS case and the real and imaginary parts of ip in 
the FER case. This allows us to write C^^^ as an effective Lagrangian of two bosons 
01, 02 in both cases: 

1 22 

4°ff = 2(^A«<^l^A«<^l + ^M02 5,.02)+^01 + ^02+^^(0252 01-0l92 02), (C.ll) 
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where 

ml = 0, ml = -Agohl, u= ^1 + gohL, (C.12) 

and 

ml = ml = A , u = CO , (C.13) 

for the BOS and FER cases, respectively. The one-loop contribution is half the 
logarithmic determinant of the linear operator corresponding to >Cgfj : 

Taking the derivative of ( P.14|) with respect to the external fields and then integrat- 
ing it again from zero, combined with dimensional analysis in n dimensions finally 
gives 

jrio) _ J f(o) _^ f dy {m\ + ml) p" + 2m\ml + Au^pl 



We have seen in Section 3 that the 1-loop result in the FER case, (|3.23| ), is a 



function of the sum hi + hji only, similarly to the classical term ( |3.2ClD . The free 
energy calculated by the TBA method also has this property. It is a very non-trivial 
fact that this behaviour can be reproduced by the perturbative calculation. In fact, 
we will now show that in the diagonal case the free energy depends only on the sum 
^L + hji to all orders of perturbation theory. 

To show this, we introduce the shifted derivatives 

d2il> = d2ip + ujip , and 821!^* = d2ip* — ujijj* . (C.16) 

Using this notation, the gauged Lagrangian can be compactly written as 
2{hL + hR)'{l + go) , .,, A—— r m'd: 



Co- ^''^"- ' '7 ^^ ' "'^^ +zMv/Ao(l + V '"' ' 



+ 



XogoM l^p {^d2r - rd2^) , M2|^|2 



12-^ 




2 (l + fl^P) (l + fl 

where Cq stands for the original Lagrangian, (|C.3|) , with all d2ip and d2ip* substituted 



by the shifted derivatives ( |C.16| ). The advantage of using the form of the Lagrangian, 
( |C.17|) , is that the explicit dependence on the external fields hi, hji occurs only 
through the combination M = {1 + go){hL + hn). In addition to this, there is 
still an implicit dependence on the external fields through u hidden in the shifted 
derivatives ( |(J.16| ). This disappears, however, since the Lagrangian ( |C17D is in fact 
equivalent to another one, denoted by Cm, differing from £ only by replacing all 
shifted derivatives, 82, by ordinary ones. 
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The fact that L and Cm are completely equivalent can be best illustrated in the 
case of the 0{\) contribution to the free energy. The relevant piece of the Lagrangian 
is (|C.10|) , which can also be written as 



£(0) = 9^^* ~d^i) + ^9^0 9^0 + W \i,f . (C.18) 

The substitution 82 -^ 82 corresponds to the shift p2 ^ P2 + i^ in momentum space. 
More precisely, when integrating over p2 in Eq. ( p.21| ) the above shift corresponds 



to deforming the integration contour parallel to the real axis. This deformation is 
allowed since the singularities of the integrand at 



P2 



iu ± i^pl + M2 , (C.19) 



lie outside the region between the real axis and the shifted contour (guaranteed by 
the inegality (|3.10|)). After shifting the integration variable Eq. (|3.21|) takes the 



simple form 

nJ (27r)"p2 + M2' ^ ' 

which is easily evaluated and yields ( |3.23| ). Now Eq. ( |C20D can be derived from the 
simplified Lagrangian 

Cfl = d,r d.ij + \d,(t> d,(t> + M^ |V;p , (C.21) 

exhibiting the one-loop equivalence between C and Cm- 

It is easy to generalize this argument to show that there is a one-to-one corre- 
spondence between the Feynman diagrams generated by C and Cm and that the 
perturbative results are identical for these two models, diagram by diagram, to all 
orders of perturbation theory. This equivalence follows from the fact that the shift 
P2 -^ P2 + i^ can be performed simultaneously in all momentum variables corre- 
sponding to the lines of a multi-loop Feynman diagram. In this way we have shown 
that all dependence on the external fields comes from the M-dependence of Cm- 
This means, in particular, that the free energy in the diagonal case takes the form 
( p.24| , |3.25| ) and it is a function of the sum h = h^ + hji only, which is exactly the 
same as the free energy for the f/(l)R case (as a function of hn). 

Appendix D. Renormalization Group Equations 

Due to its SU(2)lxU(1)r symmetry, the model (|1 . 1| ) is renormalizable in per- 



turbation theory. (Of course, in Eq. (|LT|) actually the bare coupling Aq and the bare 
f^-parameter 5^0 occur.) Using dimensional regularization, the relation between the 
bare parameters, (Ao,5'o), and the renormalized ones, {X,g), is written as: 

Xo = fi'\Zx, and I + go = {1 + g) Zg . (D.l) 
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In Eq. ( p.l| ) the dimensionful parameter, /x, is introduced as usual, to carry the 
mass dimension of Aq (which is dimensionless in 2 dimensions) . The renormahzation 
constants contain only pole terms in the variable e = 2 — n: 

Z, = l + £^. (D.2) 

Z, = l+£!^. (D.3) 

r=l ^ 

The residues of the poles, yr and Wr, can be calculated in perturbation theory: 

yr{X,9),Wr{X,9) = 0{y). (D.4) 



It is clear from ( p.l|) that ii g = —1, then also go = —1. This is a conse- 



quence of the fact that at this value of the deformation parameter the model ( |1 . 1|) 
reduces to the 0(3) nonlinear cr- model (plus a free boson), which is obvious from 
the parametrization used in Section 4. An other special value of the deformation 
coupling is g = implying go = and Wr{X,0) = in Eq. ( p.3|) . In this case the 
deformation actually vanishes and we have the SU(2)lxSU(2)r (or 0(4)) symmetric 
principal a-model. In a sense the deformed model ( |1 . 1|) interpolates between the 
0(3) and 0(4) nonlinear cr- models. 

Physical quantities depend on the renormalized couplings A, g and the dimen- 
sionful parameter /j, in such a way that the action of the renormahzation group (RG) 
operator 

o o o 

^ = /" ^ + /9a(A, g)-^ + /?,(A, g) — , (D.5) 

vanishes on them. Here the two /5-functions are obtained from the residues of the 
first order poles in ( p. 2] ) and ( p. 3D as 



/3a = A2^, and p^ = ^i + g)X^. (D.6) 



For later use we note that in the PER case the free energy, Eq. (|3.25|) , can be 
written as 

Hh) = -2h' |i±^ - ^i±^ S{go, Xo[h{l + go)r)\ , (D.7) 



where 5* stands for the infinite sum in (|3.25|) . Prom (p.7| ) we can conclude that at 
g = —1 Zx = Zg and in particular, 

u;i(A,-l) = yi(A,-l). (D.8) 
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In both BOS and FER cases, the free energy takes the form 

J=-(h) = -h"^ Fo{X,g,iJ,,h), 



(D.9) 



where the function Fq is RG invariant: T>Fq = 0. As we are interested in the 
asymptotic expansion of the free energy for large external fields we write 



h = hne* 



(D.IO) 



where Hq is a fixed, finite value of the external field and t ^ oo. Standard RG 
considerations then show that 



J^{h) = -h^Fo{X{t),g{t),f^,h 



(D.ll) 



where the running coupling A and the running deformation parameter g satisfy the 
following set of differential equations and initial conditions. 



§^^.(A,S 



dt 



MX 9) 



A(0) = A, 



(D.12) 
(D.13) 



In the deformed principal model the /^-functions are known to be of the form 

X2( ^ 



14. 15 



/3a 



47r ^ 
A^(l + 



27r 



I - g + Xv2{g) + X^vM + ■ ■ 
\l + \qi{g) + X^q2{g) + -- 



(D.14) 
(D.15) 



where the two-loop /3-function coefficients are given by |]I4|, |I5 
1 - 2^ + 5c/2 , _ 1 



P2{.g) 



and 



(lii.9) 



9 



(D.16) 



:57r ' 47r 

The three-loop coefficients P3{g), q2{g) are not known at present, but fortunately we 
will not need their explicit form beyond the combination 



M = 2g2(-1)-P3(-1). 



(D.17) 



We are interested in solutions of ( p.l2|) and ( p. 131) with asymptotically free (AF) 
behaviour, that is A ^ together with g —>■ gi as t ^ oo, where gi is some constant. 
Analyzing Eqs. ( P.14|) , (p.l5|) , we find three different AF solutions: 



1. 


91= 


9 = 90 = 9 = 0, 


(D.18) 


2. 


91 = -I 


9 = 90 = 9 = -^, 


(D.19) 


3. 


91 = -I 


-l<g<0. 


(D.20) 
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Solutions ( p.l8| ) resp. ( p.l9| ) correspond to the 0(4) resp. the 0(3) NLS model, 
while the solution ( p.20|) corresponds to the 'generic' case of deformed principal 
model. In what follows we shall concentrate on the generic case. We introduce the 
RG-invariant combination of the two couplings: 

p = 2vr hm i±fl , (D.21) 

*^oo X{t) 



which we identify with the parameter p in the S matrix ( |1.4|) . A RG-invariant 
quantity can almost be treated just as a numerical parameter. We can express 
the running deformation coupling, g, in terms of the running coupling, A and the 
RG-invariant p as 

-g = T{X,p) (D.22) 

and define an effective /3-function for A as 

Aff(A,p)=/?A(A,r(A,p)). (D.23) 

Using the perturbative expressions ( p.l4| ) and ( p.l5| ) together with ( p.l6| ) and 
( P.17|) we have 



p 
where 



T{\,p) = -1 + -f A + p(2)A2 + p(3)A3 + . . . , (D.24) 

/TT 



r)(^)--^ and r3)_ P'(P + 2) pu 

P - 8^2' and p - g^^3 +g^. (D.25J 

The unknown three- loop coefficients contribute only to u defined in Eq. ( p.l7|) . It 
is not difficult to show using Eq. (|D.8|) that 

2q2{-l) = psi-l) , (D.26) 

implying u = 0. 

Using the perturbative form ( p.24| ) in the definition ( p.23| ) we have 

PM,P) = -^ + ^A^ + . . . , (D.27) 



which allows the definition of the RG-invariant A-parameter in the MS scheme the 
usual way: 

Ams = /^^"^(^)^'^'^ ^"{l + ^(^)} • (D-28) 

The next step in the RG analysis is the introduction of the effective coupling 
A„ff(/i) defined by the transcendental equation 

|^ + (|-l),„|^^l.^. (D,29) 

^off ^ ^off ^1-MS 
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The advantage of using the effective charge is that it is a function of the physical 
quantity 

s = ln-^, (D.30) 

only and moreover the running coupling can be expressed in terms of Aoff (/i) pertur- 
batively (in the sense that it is an infinite power series): 



X = X., + ^{\n^-l)xl + --- . (D.3i; 



27r ^ fi 

The asymptotic expansion (for large s) of the effective coupling (containing terms 
oc In s) can be written as 

27r 7r(» - 2) , , , 

s s^ 



So far we have concentrated on the generic case (D.20). The 0(3) symmetric 



limit, case ( p.l9|) , {p -^ 0) corresponds to putting simply p = in Eq. (|D.32| ). The 



0(4) limit, (p.l8|) , however, cannot be obtained from the generic case although it 
corresponds formally to the p — >• cx) limit. (See the discussion in Subsection 2.3.) 

For lattice regularization some of the above formulae are modified. First of all, 
instead of ( p.l|) and (p.2| - |D.3| ) the relation between the bare parameters (Aq, ^o) and 
the renormalized ones (A, g) are given by 

Ao = AZa, and 1 + go = (1 + g) Zg , (D.33) 

where 

oo _ 

Zx = l + 5:y.(A,^)(ln/ia)^ (D.34) 

r=l 

oo 

Zg = 1 + ^ w^(A,^)(ln/ia)^ (D.35) 



r=l 



Here the parameter fl is the lattice analogue of the /x-variable of the dimensional 
scheme and naturally it is /2 that appears in the lattice version of ( p.5|) . The lattice 
/^-functions are given by 

Px = -\yiCX,g), and Pg = -{l + g)wiC\g). (D-36) 

The rest of the RG analysis is identical to the one in the dimensional scheme. 
Although the higher order /3-functions are different, the main features of the RG tra- 
jectories ( p.l8[ - [[X^UD remain the same and also the first two coefficients oi (3cs{X,p) 
are unchanged. The lattice version of the A-parameter is given by 

A. = *ie-T(l)<«-''{l + 0(A)}. (D.37) 
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